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Abstract. We propose a new lifting and recombination scheme for rational 
bivariate polynomial factorization that takes advantage of the Newton poly- 
tope geometry. We obtain a deterministic algorithm that can be seen as a 
sparse version of an algorithm of Lecerf, with now a polynomial complexity in 
the volume of the Newton polytope. We adopt a geometrical point of view, 
the main tool being derived from some algebraic osculation criterions in toric 
varieties. 



1. Introduction and main results 

This article is devoted to develop an algorithm for factoring a bivariate polyno- 
mial / over a number field K by taking advantage of the geometry of its Newton 
polytope. Geometrically, this corresponds to decomposing the curve defined by / in 
a suitable toric surface X . We will thus talk about a toric factorization algorithm. 
The usual case of dense polynomials corresponds to the classical projective comple- 
tion X = P 2 of the complex plane. Our approach is based on algebraic osculation. 
The central idea is that we can recover the decomposition of the curve C C X 
defined by / from its restriction to a suitable toric Cartier divisor D. In a previous 
work |27j . we developed a similar method based on vanishing-sums criterions and 
obtained an exponential complexity toric factorization algorithm. In contrast, we 
use here a lifting and recombination model based on a vector space basis compu- 
tation which conduces to a polynomial complexity algorithm. Our method can be 
regarded as a toric version of the algorithms developed by Lecerf [18] , [19] and by 
Cheze and Lecerf [5] for dense polynomials. Let us expose our main results. 

Main results. Let K be a number field and let / € K[ti,t2] be a bivariate 
polynomial. Suppose that / has monomial expansion 

/(*) = E c ™ fn < 

m£N 2 

where m = (mi,m 2 ) and t m = t™ 1 ^ 2 . The Newton polytope Nf of / is the 
convex hull of the exponents m for which c m is not zero. An exterior facet F of 
Nf is a one-dimensional face whose primitive inward normal vector has at least 
one negative coordinate. The associated facet polynomial of / is the univariate 
polynomial obtained from fp — X^meF c ™X m after a suitable monomial change of 
coordinates. In all of the sequel, we assume the following hypothesis 



(H\) The polytope Nf contains the elementary simplex of] 
(H2) The exterior facet polynomials of f are squarefree. 
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We denote by u) the matrix multiplication exponent. It's well known [15j that 2 < 
to < 2.37. In all of the sequel, rational factorization means irreducible factorization 
over K. Our main result is the following 

Theorem 1. There is a deterministic algorithm that, given f G K[fi,t2] which 
satisfies (Hi) and (H2), and given the rational factorization of the exterior facet 
polynomials of f , computes the rational factorization of f with C(Vol(iV/) w ) arith- 
metic operations in IK. 

In some cases, our complexity improves that of the fastest actual algorithms 
which would treat / as a dense polynomial. For instance, if Nf is the convex hull of 
{(0, 0), (a, 0), (0, b), (n, n)} for some fixed small integers a, b, we obtain a complexity 
0(n u ) while the fast algorithm of Lecerf [19] would have complexity 0(n UJ+1 ). In 
any case, the degree sum of the exterior facet polynomials of / is smaller than the 
total degree of / so that the unavoidable univariate factorization step is faster using 
the toric approach. The gain might be considerable : If a = b = 2 in the previous 
example, we use two univariate factorizations in degree 2 instead of one univariate 
factorization in degree 2n. 

Let us explain the main tools for proving Theorem 1. By the hypothesis (Hi), 
we can consider a complete regular fan X that refines the normal fan S/ of Nf 
and that contains the regular 2-dimensional cone generated by the canonical basis 
of R 2 . Such a fan determines a smooth complete toric surface X = and a 
torus-equivariant embedding of the affine plane C 2 = SpecC[ti, £2] into X. The 
rational factorization of / correspond to the decomposition over K of the Zariski 
closure C C X of the affine curve defined by /. The geometry of Nf is related to 
the intersection of C with the boundary divisor 



of the toric completion X, and we want to use this information. 

Let Div(X) be the group of Cartier divisors of X. We definitively fix D £ Div(X) 
effective with support \D\ — \dX\. For convenience, we identify D with the induced 
subscheme (\D\,Od) of X and we denote by Div(D) the group of Cartier divisors 
of D. The inclusion morphism i : D — > X induces a restriction map i* on the 
subgroup of divisors of X who intersects D properly. In particular, we can consider 
the restriction 



of C to D. The main idea is that for D chosen with sufficiently big multiplicities, 
we can recover both the rational and the absolute (over K) factorization of / from 
7c- 

The irreducible decomposition of C (~1 dX over K is indexed by the set V of the 
monic irreducible rational factors of all of the exterior facet polynomials of / and 
we decompose 7c accordingly as 



where jp corresponds to lifting P to a local factor of / modulo a local equation 
of D (see Subsection 3.1). The recombination problem consists in computing the 



dX := X \ C 2 



7c := i*{C) 6 Div(D) 
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partition of V that corresponds to the rational decomposition of C. To this aim, 
we introduce the free Z-module 

V z := { ^ eZ } c Div(X>), 

Per 

and the submodule 

V Z (D) := {7 G Vr, 3Ee Div(X), i*(E) = 7} 

of divisors of D that extend to X. We set V := V% <8> K and V(D) := Vz(D) ® K. 
By construction, the irreducible rational decomposition C = Ci U ■ ■ ■ U C s of C 
generates a vector subspace 

(71,..., 7.) CV(D) 

where 7, := i*(Cj). By the hypothesis (#2), the j/s are pairwise orthogonal in 
the basis {'lp)pev of and so dim V(_D) > s. The following theorem asserts that 
equality holds for D big enough. We say that a basis {v\, . . . , v n ) of V(D) C is a 
reduced echelon basis of V(D) if the matrix with j th row i/j is in its reduced echelon 
form in the canonical basis of the input space V (see [24] ). Such a basis exists and 
is unique. We obtain the following 

Theorem 2. Let divoo(/) be the polar divisor of the rational function of X induced 
by f . If the inequality 

D>2dW oa (f) 

holds, then (71,. .. , 7 S ) is the reduced echelon basis ofV(D). 

The proof consists in associating to 7 € Vz(D) a rational 1-form with polar 
divisor controled by C. For D big enough, that form is closed and a theorem of 
Ruppert [23] combined with a Galois theory argument permits to conclude that 7 
is Z-combination of the 7j's. There are examples in the dense case that show that 
the precision D = 2div 00 (f) in Theorem 2 is asymptotically sharp (see [19J. 

In order to apply Theorem 2 to the factorization problem, we need to determine 
an explicit system of equations that gives the vector subspace V(D) C V. To this 
aim, we use a theorem of the author that characterizes the lifting property. We 
show in [27j that there exists a morphism 

* : Div(D) <g C -> H°(X, n 2 x (D)) v 

so that 7 £ Div(D) extends to X if and only if ^(7) = 0. Roughly speaking, the 
linear form ^(7) sends a rational form u> S H°(X, Q 2 X (D)) to the sum of residues of 
a primitive of to along a local analytic lifting curve of 7. In some sense, this result 
can be regarded as a converse to the classical residue theorem, we refer to [27] for 
details. This permits to prove the following 

Theorem 3. Suppose that D satisfies the hypothesis of Theorem 2. Then, 

V{D) = ker(vl) 

for some explicit matrix A = (ap, m )pEV,mEM with coefficients in K, where M is 
the set of interior lattice points of the polytope 2Nf. 
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So we can solve the recombination problem with linear algebra over K. Then, 
we compute the rational factors of / by solving systems of affine equations. We 
finally obtain a deterministic polynomial complexity algorithm for rational toric 
factorization of bivariate polynomials. We describe briefly the main steps of the 
algorithm. The given complexities are obtained in Corollaries 1 and 2 in Section 3. 
As in [15], we use the notation O for the soft complexity. 

Toric Factorization Algorithm 

Input: / 6 K[£i,i:2] satisfying hypothesis (-Hi) and (i?2)- 
Output: The irreducible rational factors of /. 

Step 0: Univariate factorization. Compute the set V of irreducible rational factors 
of the exterior facet polynomials of /. 

Step 1: Lifting. This is the 7p's computation step. For each P € V, compute the 
associated local factor of / modulo the local equation of D = 2 div^/). This step 
has complexity d(Vo\(N f ) 2 ). 

Step 2: Recombination. 

a) Build the matrix A of Theorem 2. This step has complexity 0(Vol(Nf) 2 ). 

b) Compute the reduced echelon basis associated to A. This step has complexity 
0(Vol(7V / )Card(7') w - 1 ). 

Step 3: Factors computation. Solve some affine systems of linear equations over K 
to recover the rational factors of /. This step has complexity 0(Vol(A/) u; ). 

As already mentioned, the great advantage of our algorithm is that it replaces 
the usual univariate factorization in degree d = deg(/) by the factorization of the 
exterior facet polynomials of / : their degree sum is at most d, and much smaller in 
many significant cases. It follows that both the number of unknowns and equations 
in the recombination process decrease too and the basis computation step 2 b) is 
faster than in [TS]. Steps 1 and 2 a) rely on classical modular algorithms (Newton 
iteration, modular multiplication) whose complexity analysis is delicate due to the 
sparseness of /. This partially explains that our lifting complexity does not reach 
the soft complexity 0(d 2 ) obtained in [19] for dense polynomials. Step 3 has the 
highest cost of the algorithm because in the general sparse case we have to use 
linear algebra instead of the fast multiplication or partial fraction decomposition 
methods that are used for dense polynomials ([H], [H], [H])- Finally, let us mention 
that the algorithm developed by Lecerf [18] in the bidegree case suggests that it is 
possible to reduce both the number of facet factorizations and the lifting precision. 
We refer to Subsection 3.5 and Section 4 for further comments. 

Related results. Classical results about polynomial factorization can be found 
in |15j . For more recent advances, we refer the reader to the introduction of [8] 
(and to the complete list of references therein) that gives a large and comprehensive 
overview of the current algorithms for factorization of polynomials. We only discuss 
here the most related results. 

Using linear algebra. Factoring multivariate polynomials by means of linear alge- 
bra has been made possible by the powerful irreducibility criterion of Ruppert [23j . 
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This is the so-called logarithmic derivative method, that relates the basis compu- 
tation of the vector space of closed rational 1-forms with some appropriate polar 
divisor. This point of view has been developed by Gao in [14j . who combined the 
logarithmic derivative method with the Rothstein-Trager algorithm for absolute 
partial fraction decomposition ([15] . Theorem 2.8). Finally, as pointed out in the 
introduction, Lecerf [19], [18] and Cheze-Lecerf [8] recently developed very efficient 
hybrid algorithms for rational and absolute factorization, by combining Gao's ap- 
proach with a lifting and recombination scheme. This is the point of view we follow 
here. 

Using Newton polytopes. Factoring polynomials by taking into account the Newton 
polytope is an active area of research. In llj, M. Elkadi, A. Galligo and the author 
use some probabilistic interpolation criterions [26] , by replacing the divisor D with 
a generic ample curve "close to the boundary". In [57], the author looks for the 
effective decompositions of jc that may be lifted to X. By taking into account 
natural degree conditions imposed by the Minkowski-sums decompositions of Nf, 
there appears supplementary vanishing cohomology properties of the osculating 
divisors that permit to use the smaller precision D = div oc (/) + dX (see also 
[18] for a similar comparison in the dense case). In return, it gives a problem of 
partitions of V(D) D {0, 1} V that has exponential complexity in the worst case (see 
Subsection 3.6). A comparable algorithm is obtained in [I], where the authors use a 
more combinatorial approach. In [3], the authors show that the low degree factors 
of / can be computed in polynomial time with respect to the fcwnomial encoding 
of /. 

Organization. We prove Theorem 2 in the next Section 2. In Section 3, we develop 
a toric factorization algorithm and we prove Theorem 1. In Section 4, we compare 
our method with the most related dense and toric algorithms and we discuss some 
possible improvements. We conclude in Section 5. 

2. Proof of Theorem 2 

We follow the notations of the introduction. We saw that the rational decompo- 
sition C = CiU- ■ -UCs of C generates a vector subspace (71, ... , 7 S ) C V(D) and we 
want to show that the opposite inclusion holds when D > 2diVoo(/)- The strategy 
consists in associating to 7 € Vz(-D) a closed rational 1-form to on X whose polar 
divisor is controled by C. A theorem of Ruppert [23 implies that to is a C-linear 
combination of the logarithmic derivatives of the absolute factors of /. Finally, we 
conclude by Galois theory that 7 is Z-combination of the 7j's. 

We need first two preliminaries lemmas that clear up the behaviour of restriction 
with respect to derivation. The remaining part of the proof will follow in Subsection 
2.2. If not specified, all schemes are considered over C. 

2.1. Notations and preliminaries lemmas. We denote by Id the ideal sheaf of 
D and by Od its structural sheaf. The structural sequence of D is 

(1) -> I D -» Ox £ O d -> 0, 

where the restriction map i* is induced by the inclusion i : D — > X. We denote by 
Ox{D) the sheaf of rational functions with polar divisor bounded by D, by il q x the 
sheaf of regular g-forms and we let Q q x (D) := Q x ® Ox{D)- 
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We say that B G Div(X) is a normal crossing divisor if it has local equation 
x\ ■ ■ ■ x r = where the Xi's form part of a local system of coordinates (xi, . . . ,x n ) 
of X (so n = 2 in our case). For such a -B, we introduce the sheaf n q x (logB) of 
rational g-forms with logarithmic poles along B. By definition, (j> G f2^(log-B) if 
and only if both h<f> and hdcp are regular for some local equation h — of B. It is 
well known that fl q x (logB) is a locally free sheaf of Ox-module [25j . 

The following lemma clears up the behaviour of the restriction morphism with 
derivation. 

Lemma 1. Let B, D as before, with \D\ C \B\. Let F be an effective divisor which 
intersects D properly. The differential d induces a commutative diagram 

tl x (\ogB)®O x (F) A n\{B)®O x {2F) 

li* [i* 
VL x {\ogB)®0 D {F) H il x {B)®0 D {2F). 

Proof. We show Lemma 1 for an arbitrary smooth complete variety X of dimension 
n. Since B is normal crossing, it has local equation x\ ■ ■ -x r = where the x^a 
form part of a local system of coordinates (x\, . . . , x n ) of X. The sheaf ft x (logB) 
is a locally free sheaf of Ox-modules and a germ <p G Q x i}°gB) <8> Ox(F) has a 
unique representation 

4> = h\dx\jx\ + • • • + h r dx r /x r + h r+ idx r+ \ + • • • + h n dx n 

for some hi G Ox{F) (see [25], p. 186). It is clear that dhi has its polar divisor 
bounded by 2F. We deduce that 

d<f) = dh\ A dxi/xi + ••• + ••• + dh n A dx n 

belongs to Q X (B + 2F) and the upper row is well-defined. In order to show that 
do is well-defined, we need to show that 

4> G n x (B + F) <g) I D =>d<j)£ n x (B + 2F)®I D . 

Since D is supported on \B\, it has local equation x k := x kl ■ ■ ■ x kr for some k, L G N. 
Thus if <fr G £l x {B + F) ® I D , we have hi = x k ti i for some h[ G O x {F) and 

dhi dxi , . \ -> , dxj dxi „ , dxi 
—± A — - = h, kj — - A — - + dh\ A — - 

3=1 J 

belongs to Q X (B + 2F) for all i = 1, . . . , r. In the same way, it is easy to check 
that dhi A dx l /x k G Q X (B + 2F) for i > r. Multiplying by x k , we obtain that 
d(f> G n 2 x (B + 2F) ® I D - □ 

We now pay attention to the behaviour of the restriction map with logarithmic 
derivation. We denote by M.x,d the sheaf of rational functions whose polar locus 
intersects D properly. We have an exact sequence 

(2) -> I d M x ,d -> Mx, D ^M D ^0, 

where Md '■= M.x,d ® Od is the sheaf of rational sections of Od- The multiplica- 
tive version of (2) is 

(3) 0^ 1 + I d M x ,d ^M* xd ^M* d ->0, 
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where * stands for the multiplicative sheaves of units. On the other hand, the 
logarithmic derivative dlog(h) := dh/h induces the natural morphisms 

(4) (flog : M* XtD -^n x ® M x ,d and d\og : 0* x -► Q x 

of sheaves of abelian groups. We have the following 

Lemma 2. Let B, D as before and suppose that \D\ C \B\. The morphisms in (4) 
combined with the natural inclusion j : Q x ® Mx,d —> fix (log B) ® .A4x,d induce 
the commutative diagram 

M* X>D £> Mb - 

Ijorflog jcfolog 

n x (iogB)®M x , D n x (io g B)®M D -► o, 

and its regular version 

o* x ^ o* D -► 

I I 
VL x {\ogB) ^ n x (\ogB)(g>0 D 0. 

Proof. Let a e .MJ, and u 6 ,M X n so that a = i*(u). By (4), the morphism 

d D log : M* D £l x (\ogB) Mb, cfo log(a) := o dlog(u)) 

will be well-defined if we show that 

u G 1 + IdM x . d => d D log(u) = 0. 

So let it = 1 + ft, for ft e Id-Mx.d a germ at some smooth point of £>. Thus B and 
_D have respective local equation x — and x fc = for some k > 0, and ft = mi' 



i v vj^i in i la' i i .( — ■■ ■ nuu .< — \| iv'i " ■ i i n ■ ■> < 1 1 j.v i / -■ — 1 1 tX 

where m E M.x,d- So 

„ , . kmx k ~ 1 dx + x k dm mx k ,dx dm, 

dlog 1 + ft) = — -j = - fc— + 

1 + xm 1 + a; K m a; m 

belongs to the subsheaf Q x (\og(B)) <g> LdM x ,d C fi x <8> A^x,r>- By tensoring (2) 
with the locally free sheaf fl x (\og(B)) we deduce that 

i*(j o dlog(l + ft)) = g O x (logB) ®M C . 

The divisor B being normal crossing, we check easily that the same conclusion 
holds when ft is a germ at some singular point of B. This implies that the map 
dp log is well-defined, giving the first diagram. The regular version follows by 
letting m E O x in the previous reasonning and by using the multiplicative version 
of(l). □ 



2.2. Proof of Theorem 2. We come now to the proof of Theorem 2. We denote 

by 

T := SpccC^f 1 ,^ 1 ] and C 2 := SpecC[ii, t 2 ] 
the complex torus of X and the affine plane endowed with canonical coordinates 
t = (ti 7 t 2 ). We identify rational forms of T and C 2 with the rational form they 
induce on X. We suppose from now that D = 2 div<x>(/) and that B = X \ T. The 
toric surface X being smooth, the toric divisor B is normal crossing. Since / has 
no poles in the torus, we have \D\ C \B\. 
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Let 7 G V(D). We need to show that 7 is linear combinations of the 7/s. There 
is no loss of generality to suppose that 7 G Vz(D). Let T := C (~1 <9X be the 
intersection of C with the boundary. By hypothesis, we know that 

(5) 1 = ^ jJlplp and 7 = i*(C 7 ), 

per 

where 7 P G Div(D) is induced by the germ of C at p, [i p is an integer and where 
C y G Div(X). Since the torus has a trivial Chow goup, there exists E 1 supported 
on B so that 

div(g) = C 7 - E-y 

for some rational function g G C(X). The following key lemma ensures that the 
poles of the restriction i*(dg/g) are controled by C. 

Lemma 3. We have i*(dg/g) G (log 5) ® O d (C)). 

Proof. Obviously, i*(dg/g) defines a rational section of Sl x (logB) ® C?_d and we 
need to show that the germ g p of g at p satisfies 

i*(dg p /g p ) £ n x (log B)®O d (C) 

for all p € \D\ (for convenience, we omit the index p in the stalk notations). Let us 
write g p = G p /H p for some local equations G p and H p of respectively C 7 and £' 7 . 
Thus 

(6) dg p /g p = dG p /G p - dH p /H p . 

Since E 7 is supported on B, we have dH p /H p G (logS) and it's enough to show 
that 

i*(dG p /G p ) e O^(logB) ® D (C) 

for all p G |D|. For convenience, we let 7 P := and [i p := for p G |D| \ |T|. By 
(5), the germ (C 7 ,p) of C 7 at p satisfies 

i*(C 7 ,p) = fi p j p = n P i*(C,p) 

for all p G \D\. This is equivalent to that 

(7) i*{G p /F^)eO* D , 

where F p is any local equation of C at p. Lemma 2 combined with (7) implies that 

(8) d D log(i*(G p /Fg»)) = i*{dG p /G p ) - v P i*(dF p /F p ) e ^(log-B) ® O d . 

Since i*(dF p /F p ) G ^(logB) <g> O d (C), it follows that i*(dG p /G p ) belongs to 
^(logS)®C D (C). □ 

Lemma 4. T/iere errasfc a unique rational form u G H°(X, 0^-(log(B) <g> 0x(C)) 
suc/i i/iat i*(dg/g) — 

Proof. By tensoring (2) with the locally free sheaf fi^- (log B) ® Ox (C) and by look- 
ing at the associated long exact cohomological sequence, we deduce from Lemma 3 
that it is enough to show that 

(9) H 1 (X ,Q x (log B)®O x (C)® 1 D ) = 0. 
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But C — divoo(/) = div(/) being principal, multiplication by f 2 gives a global 
isomorphism 

O x {C)®I D ~O x {-C). 
Moreover, we know by [T3] that the sheaf H, x (\ogB) is globally trivial. Thus, there 
is a global isomorphism 

n x (logB) <g> O x (C) ®I D ~ Ox(-C) © O x {~C). 

Since C has a non negative intersection with each irreducible toric divisor of X, 
the line bundle Ox(C) is numerically effective [13]. Moreover, it is well known [13] 
that we have equality 

(10) H°(X, O x (C)) = {t m /f, m e N f n Z 2 }. 

By (Hi), the Newton polytope Nf has dimension 2, and we deduce from (10) that 
Ox{C) is big. It is a standard result that 

H\X, O x (-C)) = 

for any big and nef Cartier divisor C on a smooth complete variety X (see [17] . 
Theorem 4.5 for instance). Finally, (9) holds and there exists u> as in Lemma 
4. Since H°(X, Ox{— C)) — 0, the previous reasonning implies that such an uj is 
unique. □ 

Up to here, we can show that all previous lemmas would remain valid with the 
choice D = diVoo(/) + dX used in [57] • The choice D > 2div oc (f) appears to be 
essential in order to have 

Lemma 5. We have dui = 0. 

Proof. Lemma 1 applied with F = C combined with Lemma 4 gives 

t*(dcj) = d D (i*(uj)) = d D (i*(dg/g)) = i*(d(dg/g)) = 0, 

where i*(dw) G H (X,fl 2 x (B + 2C)®O d ). By tensoring (2) with Q 2 X (B + 2C), 
and by using the associated long exact cohomological sequence, we deduce that 

duj e H°(X, fl x {B + 2C) ® I D ). 

Since 2C - D = div(/ 2 ), it follows that f 2 duj E H°(X, U X (B)). By [T^, the divisor 
B = X \ T is an anticanonical divisor of X and there is an identification 

H°(X,n x (B))=' 



, dt\ A dt 2 



Thus, 

, c dti A &2 

00) = — 

J 2 ht 2 

for some constant c G C. This form being exact, its residue at zero vanishes. This 
forces c// 2 (0) = res {diu) = (recall that /(0) ^ by (H t )) and = 0. □ 

Here comes a theorem of Ruppert in the picture. 
Lemma 6. There exists some constants Cj,a\,a2 S C such that 

Edq~j dti dt2 
cj— + ffli- h a 2 — , 
. , ffj ii *2 

with , <?t t/ie irreducible absolute factors of f. 
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Proof. By [13], the map in i— » dt m /t m gives an isomorphism Z 2 (gi Ox — fix (log £?). 
We deduce that the map 

H°(X,O x (C))®H°(X,O x (C)) -» iT°(X,fi x (lo g B)®0x(C)) 
(ri,r 2 ) h-> ndti/tx + ridti/ti 

is an isomorphism. It follows from (10) that there exists (unique) polynomials 
hi , h 2 such that 

(11) ^Tr + TT' ^CJVy. 

J '1 J 12 

On the other hand, cj being closed by Lemma 5, its restriction to C 2 defines an 
element of the first algebraic De Rham cohomology group iJ x (C 2 \ Co), where 

Co := (C + B)nC 2 = {t 1 t 2 f = 0}. 

By a theorem of Ruppert |23j , it follows that there are uniquely determined con- 
stants ci, . . . , Ct, oi, ai € C and a unique ezaci rational 1-form u>' such that 

, r—s d(Ji dt\ dt 2 

ui — uj = y Cj^-+a\ h a-2 — • 

f-f ?j *i *2 

2 — 1 J 

Since the right hand side can be written as in (11), we deduce that there are 
polynomials pi,P2 so that 

, pxdti +p 2 dt 2 , , , . 

w = — — ; , deg(^) < deg(/*it 2 ), 

1112/ 

where deg(-) stands for the total degree. Since the polynomial t\t 2 f is reduced (by 
(Hi) and (H 2 )) and u>' is exact, it follows from [Sj, Proposition 3 that ui' = 0. This 
gives the desired expression for ui. □ 

Lemma 7. Lei C,- £ Div(X) be the component of C defined by the absolute factor 
Qj . We have the relation 7 = X^=i c i**(Cj)- 

Proof. Let p £ \T\ and let x — and y — be some respective local equations of C 
and _B at p. By Lemma 6, the germ ui p of w at p satisfies us p — Cj dy/y 6 Q x (log _B) , 
where Cj is the unique component of C passing throw p (unicity of Cj comes from 
(#2))- It follows that 

t*(wp) - Cj i*(dy/y) E n x (\ogB) ® D , 
while (6) combined with (8) implies that 

i*(d 9p /g p ) - ft,? (dy/y) G ^ x (logS) <g> C D . 
Using equality i* (dg p / g p ) — i*(lo p ) induced by Lemma 4, we deduce that 
(fi p - Cj )i*(dy/y) e n x (\ogB) ® £> D , 

so that 

(pp - Cj)i*(dx A dy/xy) <E £l x (B) (g> O D . 
Hence, there exists ip € £l x (B) such that 

(fip - Cj)dx A dy/xy -i^e^fB)® O x (C) <g> Id. 
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Since | divoo(/)| = \dX\ (see Subsection 3.1), we deduce that (D,p) > (B,p) for any 
p E |r|. It follows that the previous germ of 2-form has its polar divisor bounded 
by the smooth germ of curve {C,p) = {y = 0}. Hence it has no residue at p 

reSp [(/ip — Cj)dx A dy/xy — ip] = 0. 

In the same way, tjj e Vl 2 x {B) forces ies p (ip) — so that 

= reSp — cj)dx A dy/xy] = /i p — Cj. 

The relation 7 = Y^j=i c j**(Cj) follows. □ 

Lemma 8. 7 is Z- combination of the jj 's. 

Proof. Let Cj and Ck be conjugate components. We need to show that Cj = c^. 
Let us consider the schemes X, D , Cj and 7 as schemes over Spec IK (we keep the 
same notations for simplicity). Since both D and 7 are defined over K, the group 
Autx(^) of K-automorphisms of X acts on Div(D) and fix 7. Let a G AutR(X) 
be such that o-(Cj) = Ck- Then 

(12) c k i%C k )+Y,c i i%C i )= 1 = a( 1 ) = Cji*(C k )+J2c i o-(i*(C i )). 

Since o induces a permutation of the irreducible absolute components of the rational 
curve C, (12) implies that (cj ~ Ck)i*(Ck) is supported on X^fc c ii*{Ci)- Since by 
the hypothesis (H2), the schemes i*{Ck) and X^ fc «*(Ci) have disjoint support, 
this forces equality Cj = Cfe. □ 

This shows that V(D) = (71, . . . ,7 S ). Since the Tj's have coordinates in {0, 1} 
and are pairwise orthogonal in the canonical basis {"fp)pev of the ambient space 
V , they form (under some unique permutation) the reduced echelon basis of V(D). 
Obviously, this remains true for any choice D > 2divoo(/). Theorem 2 is proved. 
□ 

Remark 1. If we rather consider the 'L-module Vz induced by the irreducible decom- 
position of^/c = CnD over the algebraic closure K o/K, we can check that Theorem 
2 remains valid. Namely, if D > 2div oc (,f), then the submodule Vz(D) C Vz 0/ 
divisors that extend to X is free generated by the restrictions of the irreducible ab- 
solute components of C . This might be usefull to compute the absolute factorization 
off- 

3. A TORIC FACTORIZATION ALGORITHM 

We come now to the proof of Theorem 1. We introduce notations and review 
basic facts about toric geometry in Subsection 3.1. We describe the echelon basis 
computation in Subsection 3.2 and the factors computation in Subsection 3.3. We 
develop a toric factorization algorithm and prove Theorem 1 in Subsection 3.4. 
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3.1. Notations and preliminaries. We refer to [10] and [13] for an introduction 
to toric geometry. As before, X designs the smooth toric surface associated to a 
regular fan £ that refines the normal fan of Nf . 

We denote by Do, . . . , A+i the irreducible toric divisors of X and by po, . . . , p r +i 
the corresponding rays of S. Since £ is regular, we can order the A's in such a 
way that the generators rji of the monoids pi n Z 2 satisfy det(?7i, ?7i+i) = 1, with the 
convention r/ r+ 2 = ?/o- We denote by Ui ~ C 2 the affine toric chart associated to 
the two-dimensional cone piM. + © pi + iM. + . Thus, 

Ui = Spec C[x, y] 

where affine coordinates and torus coordinates t = (ii, £2) ar e related by relations 

for all raeZ 2 , with (•, •) the standard scalar product (see [10] for instance). More- 
over, the irreducible toric divisors have affine equations 

A n Ui = {x = 0} and A+i D Ui = {y = 0}, 

and I Dj-I n E/j = for j <£ {i, i + 1}. 

Since the fan S contains the cone generated by the canonical basis (ei, e-i) of M 2 , 
we can chose the indexation such that (r] r +i, Tjo) — (ei, 62), in which case 

dX = D x -\ \-D r . 

It is well known that the toric divisor A appears in the principal divisor div(/) 
with multiplicity 

di := -min meNf (m,rji). 
By (Hi), we have do = d r +\ = 0, and di > otherwise, so that 

divoo(/) = diA 1- d r D r 

for some di, . . . , dr > 1. In particular, | div oc (/)| = \dX\ as asserted in the proof of 
Lemma 7. Since / = J2m£N f nz 2 c mt m , the rational function 

(13) fi{x,y):= c m x< m ^+V m * +1>+<il+1 

is a polynomial which does not vanish at (0, 0), and C has local equation 

CnU l = {f t (x 1 y)=Q} 

in the chart A We call the univariate polynomial Pj(y) := fi(0,y) the i^-facet 
polynomial of /, or exterior facet polynomial of / when i ^ 0, r+1. The polynomial 
Pi has degree 

k := deg(P l ) = Card(^j.' n I?) - 1, 

the lattice length of the i*' 1 exterior face 

Nf := {m G iV>, (m,»/i) - 

of iV/. If pi belongs to the original normal fan S/ of iV/, this face has dimension 
one. Otherwise, it is a vertex of Nr. 

A point p £ C fl A has local coordinates (0, y p ) in A where j/ p G K is a root of 
Pi. Under the hypothesis (P2), the curve C intersects transversally the boundary 
of X so that there exists a unique series (j) p G such that <^ P (0) = y p and 

<j) p ) = 0. The restriction j p of the germ of C at p to an effective toric divisor 
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D = ^2 kiDi is thus uniquely determined by the truncation at order fcj of the series 
4> p . Since the <f) p 's are conjugate when the j/ p 's run over the roots of an irreducible 
rational factor P of Pi, the Cartier divisor of D 

ip ■■= E > 

P(y P )=o 

is defined and irreducible over K. It follows that the restriction 7c of C to D admits 
the irreducible rational decomposition 

~tc = E 7p ' 

Pev 

where V := V\ U • ■ ■ U V r is the union of the sets Vi of the non constant monic 
irreducible rational factors of the i th exterior facet polynomial of /. 

3.2. Computing the reduced echelon basis. Proof of Theorem 2. We now 

give a way to compute the reduced echelon basis ofV(D) C V when D = 2 div oc (/). 
So we need criterions for lifting Cartier divisors from D to X. In [27], the author 
obtain such conditions that involves the algebraic coefficients of the series cp p 's. 
From an effective point of view, we rather follow [H] and work over the residue field 

K P := K[y]/(P(y)) 

associated to each PeP. Let yp be the residual class of y. A series B G Ip[[s]] 
can be uniquely written 

B = E hVp 

where lp :— deg(P) and bj := cocS(B,y p ) G K[[x]]. For convenience, we denote by 
cocfffc(B, y p ) G K the coefficient of x k in bj. 

Suppose that P G Vi- Under the hypothesis (H2), there exists a unique power 
series (fip G Kp[[i]] such that fi(x, (fip) = and (pp(0) — yp- So <j)p is the conjugate 
class of the series <j> p associated to the roots of P. Note that <pp is invertible. For 
all integer k 7^ 0, we define 

± 
k 



and we define B°(0p) := log(0p) to be the unique primitive of <f>' p /(j>p which 
vanishes at . For all m G 1? . we define 



apn 



lp-x 

E 

3=0 



TV(P) coeff^^B^^ {<j>p),y P ), 



where i is chosen so that P G Vi and where IV (P) designs the sum of the j th power 
of the roots of P. So ap m G K. Theorem 3 follows from the following 



Proposition 1. Let M denote the set of interior lattice points of 2Nf and let A 

denote the matrix (ap m )pev,meM ■ There is equality V{D) — ker(^4). 
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Proof. By the algebraic osculation Theorem 1 in [27], we know that there exists a 
pairing 

(•, -) D : Div(-D) ® C x ft|(L>)) -» C 

such that 7 G Div(D) extends to X if and only if (7, -)rj = 0- Since D — 2 divoo(/), 
we have 



H°(X,n 2 x (D)) = C^ m , Vm := t> 



mem tlh 
and the explicit formula of Proposition 1 in |2T] gives equality 

{IpMd = coeff(B< TO '*+ 1 >(^),s-<™''»>) 

for all p G C n Z?i. Note that — (m, r?i) < 2di by hypothesis, so that the previous 
expression only depends on <\> p modulo (x 2di ). 

Suppose that y p is a root of P G V%- Since 4>p is the conjugate class of <p p , there 
exists for all k G Z a unique polynomial R fc G K[a;][y] with degree < 2di in a: and 
degree < lp in y such that 

B fe (^ P ) = R k (y P ) mod (z 2d *) and B fe (0 p ) = R k (y p ) mod (z 2 *)- 

Hence, we have congruence relations 

Ip-i ip-i 
J2 B fc (0 p ) = E coeff(R fe ,^)i/^ J] Tr^(P)cocff(B fe (0p),^) 

modulo (a; 2di ) and the relation 

(lP,^rn)D = ^ (~(p,4>m)D = aPm 
P(Vp)=0 

follows. So 7 G Vz extends to a Cartier divisor on X if and only if -fA = 0. It 
follows that V(D) = ker(A). □ 

Let us look at the algorithmic complexity underlying Proposition 1. Following 
[T5] . we use notation O for soft O, and we let 2 < u> < 2.34 be the matrix multipli- 
cation exponent. We recall that the complexity for multiplying two polynomials of 
degree d belongs to 0{d) (Schonhage and Strassen algorithm, [15]). 

Corollary 1. Suppose given the rational factorization of the exterior facet polyno- 
mials of f . We can build the matrix A with 0(Yol(N f) 2 ) arithmetic operations in IK 
and then compute the reduced echelon basis ofV(D) with 0(Yol(N f) Card^)" -1 ) 
arithmetic operations in K. 



Proof. We divide the proof in three steps. 

Step 1. Computing the <f>p 's. We use the classical Newton iteration algorithm [15] . 
In order to estimate the cost in our toric setting, we need the following 

Lemma 9. Let P G V L , 4> G Kp[[i]] and let fi G K[:c,y] as defined in (13). For 
any k G N, we can evaluate fi(x,(j>) G Kj>[[i]] modulo (x ) with 0(kVo\(Nf)) 
arithmetic operations inKp. 
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Proof. Since /j is a sum of 0(Vo\(Nf)) monomials, we can evaluate fi(x,-) at <fi 
modulo (x k ) by evaluating each of the involved monomials with 0(Vo\(Nf) log(rti)) 
operations in Kp[[x]]/ (x k ), where rii is the total degree of fi in y. All what we need 
to show is that rii is not "too big". By (13), we have 

m := deg (ft) = max (m,r] l+ i) - min (m,r] i+ \), 

" mGNf m£Nf 

Let m £ Nf. By (Hi), the polytope Conv{(0, 0), (1, 0), (0, 1), m} is contained in 
Nf. Since it has euclidean volume (mi + m 2 )/2, we deduce ||m|| £ 0(Vol(Nf)) 
for all m £ Nf, where ||m|| := \m\\ + | rrz 2 1 - There remains to estimate ||?7j + i||. 
As before, we check easily that H77H £ 0(Vol(Nf)) for all inward primitive normal 
vectors of the one dimensional faces of Nf . Let j > i be the first index for which 
the j'^-exterior face of Nf has dimension one. So we can write = arji + br/j for 
some positive rational numbers a, b. We have relations 

&det(r/j, rjj) = det(r/j, r/ i+ i) = 1 and a det(?7j, i]j) = det(?7i + i, rjj) < det(r/i, r)j), 

last inequality using that £ is regular and refines £/ (see [9]). So a, b < 1. Since 
WmW £ 0(Yol(N f )) and \\rjjW £ 0(Vol(N f )), it follows that H^+iH £ 0(Vo\(N f )). 
Since \(m,TH+i)\ £ 0(||m|| ||?7 i+ i||), we deduce n t £ C(Vol(^/) 2 ). 

Hence, the cost for the evaluation step is C(Vol(iV/) log(Vol(iV/)) operations in 
K P [[x]]/(x k ), or 0(kVo\(N f )) operations in K P (Corollary 9.7, [T5]). □ 

By using the fast modular Newton iteration Algorithm 2, [8] and by replacing 
the given evaluation cost by that induced by Lemma 9, we deduce that we can 
compute the series cf>p with precision 2di with (D(dt Yo\(Nf)) operations in Kp. 
Each operation in Kp takes 0(lp) operations in K and we have J2pev = h- 
Since 2C and D have the same Picard class, we deduce 

r r 

2dih = 2d * de &(° ■ D i) = de 8'( C ■ D) = 2 deg(C • C) = 4 Vol(JV», 

i=l i=l 

the last equality using basic toric intersection theory (see [13] for instance). It 
follows that computing all the </>p's up to the precision imposed by D has complexity 
<5(Vol(A7) 2 ). 

Step 2. Building the matrix A. Let P £ Vi- Computing the series B^ m,I ' i+1 ^ (<pp) 
with precision 2di for all m £ mt(2Nf) n Z 2 requires at most one inversion in 
K P [[x]]/(x 2d *) and 0(Vol(N f )) evaluations of monomials in K P [[x}}/ (x 2d *)[y] of 
degrees bounded by ni = 0(Vol(/) 2 ) (and possibly the computation of a primitive 
of <fi'p/(j)p modulo (x 2di )). Each operation takes O(lpdi) operations in K, giving 
a total number of 0(Vo\(Nf)lpdi) operations in K. Summing up over all P £ V, 
we obtain a total number of 0(Vo\(Nf) 2 )) operations in K for computing all the 
B fe (^p) involved in the definition of A. Then building the matrix A has a negligeable 
cost. 

Step 3. Computing the reduced echelon basis ofV(D). Since V(D) is determined by 
C(Vol(iV/) equations and Card(7 3 ) unknowns, we can compute its reduced echelon 
basis with 0(Yol(N f ) CardCP)^- 1 ) operations in K ([24 , Theorem 2.10). □ 
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3.3. Factors computation. We want now to compute the rational factors of /. 
In all of this section, we fix 7 in the reduced echelon basis of V(D) and we denote 
by q the corresponding factor /. We first compute the Newton polytope of q. By a 
Theorem of Ostrovski [22], N q is a Minkovski summand of Nf so that it is enough 
to compute the integers 

e, := -min m&Ng {m,r]i}, i = 0,...,r + 1. 

Suppose that 7 = J2p£-p HpJp- We define 

U(Tf) ■= E ^pdeg(P) 
Per, 

for alH = 1, . . . , r. We obtain the following 

Proposition 2. We have eg = e r +i = and the integers e\,...,e r are the unique 
solutions of the affine system 

r 

^e<aij j = l,...,r, 

i=l 

with a^i+i := 1, an := det(r?j_i, 7?j+i) and ay := for j 



Proof. Since N q is a Minkowski summand of Nf, the hypothesis (Hi) forces G N qi 
so that eo = e r +i = 0, and e, > otherwise. By Subsection 3.1, we know that 

divoo(g) = e 1 D 1 + • • • + e r D r 

while the component C = div (q) of C satisfies 

(14) deg(C ■D j )=deg( llDj )=l j ( 1 ), j = l,...,r. 

Now, the Chow group of the smooth affine plane completion X is Z-free generated 
by the classes of D\, . . . , D r and numerical and rational equivalence coincide on a 
smooth toric variety [13] . By duality, it follows that (14) uniquely determines the 
class of C , that is the class of divoo(g), that is the integers e%, . . . , e r . We have 
equality 

T" 

deg(C" • Dj) = deg(div 00 (< z ) • D 3 ) = ^ e 4 deg(A • Dj) 

i=l 

and it's well known that deg(Z?i • Dj) — aij (see [13] for instance). Proposition 2 
follows. Note that det(ay) = ±1. □ 

We can now consider q(t) = J2 m eN q ni. 2 c rnt m as a vector indexed by the lattice 
points of N q . For all m £ N q and all P G V, we define 

r Pm e K P [[x]]/(x^ +1 ), r Pm := x M+e i(j) (^m +1 ) mod ^+1^ 

where i is chosen so that P G Vi. This definition makes sense since 4>p is invertible 
and (m, 7jj) + > for all m G N q . We obtain the following 
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Proposition 3. Under normalization q(0) = 1, the coefficients vector (c m ) m eN q nz 2 
of q is the unique rational solution of 

c = l and Vp m =0eKp[[i]]/(i £,+1 ) VP G V, ^ P + 0. 

mGN q nZ 2 

The induced system (S 7 ) over K contains 0(Vol(N q j) affine equations. 

Proof. Let h G K[ii,t 2 ] with Newton polytope contained in N q . So h and q define 
global sections of the line bundle Ox(E), where E := X)[=i e iA- ln particular 
h = Xq for A G C if and only if the two Cartier divisors H := div(h) — E and 
C := div(g) — E are equal. Since the restriction 

H°(X,Ox(E))^H°(X, E+dx (E)) 

is injective and H and C" are rationally equivalent to E, we deduce that 

C = H 4=» C , D(E + dX) = Hn(E + dX) 

<^=> C , D(E + dX) G H, 

with intersection and inclusion taken scheme theoretically. Since N q G Nf, we have 
E < diVoo(/), so that E + dX < D and 

C C)(E + dX) =jn{E + dX) = (J lP n(E + dX). 

Let hi = be the local equation of H in the chart f7j. For P G Pj, we know that 
7p n (P + <9X) = 7p n (ej + 1)A is contained in f7j. Hence, we have equivalence 

7p f~l (ej + 1)A C i? hi{x,<t> P {x)) = mod(a; e,+1 ). 

If = Y< m eN q ni 2 u mt m , we can chose (following (13)) 

meN q ni 2 

Since (ftp is invertible, we obtain by linearity that r )pC\(E + dX) C H is equivalent 
to that equality ^2 m u m rp m = ® holds in Kp[[i]]/(/ ,+1 ). It follows that q is the 
unique solution of the announced system of equations. 

A linear equation in Kp[[s]]/(i ei+1 ) being equivalent to a system of lp(ei + 1) 
equations over K, the total number of equations of the induced system (5 7 ) over K 
is 

r r 

1 + ^ ^ (e, + 1) Mp Zp = 1 + ^( ei + 1) deg(C • A) 
i=i p e -Pi i=i 

= deg(C" • £) + deg(C" • dX) + 1 G 0(Vol(JV,)). 

□ 

We deduce the following 

Corollary 2. Suppose given the reduced echelon basis ofV(D). We can compute all 
the irreducible rational factors of f with at most 0(Vo\(Nf) u ') arithmetic operations 
in K. 
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Proof. Let 7 belongs to the reduced echelon basis of V(D). Once the affine system 
(5 7 ) is built, Proposition 3 allows us to compute the corresponding factor q of / by 
solving a system of 0(Vo\(N q )) affine equations over K and 0(Vo\(N q )) unknowns. 
This requires 0(Vol(N q ) u ) arithmetic operations in K (Theorem 2.10 in [24]). There 
remains to build the system. Computing N q using Proposition 2 has a negligeable 
cost. Let P £ Vi- We need to compute (f>p with precision x ei+1 . Since e» < di and 
di > 0, we have + 1 < 2d^ so that <pp has already been computed with a sufficient 
precision. In the same way, computing the involved powers m £ N q from 

the already computed powers </>^ n ' 1,,+1 \ m £ int(2Nf n Z 2 ) has a negligeable cost 
too. If / = qi ■ ■ ■ q g is the rational factorization of /, we have inequality 

Vol(N gi ) + ••• + Vol(iV, s ) < VolCN,, + • • • + N gs ) = Vol(N f ) 

and we finally need at most 

0(Vol(iV 9i r + • • ■ + Vol(iV g J w ) c 0(Vol(JV/) w ) 

arithmetic operations in K for computing the rational factorization of / from the 
reduced echelon basis. □ 

Let us remark that the system (S 7 ) has a particular sparse structure. For instance, 
it contains the subsystems of type Vandermond 

£ c m y { r u+l) = 0, PEVi^P^O 

meN^nz 2 

that determines (up to multiplication by some constant) the i th exterior facet poly- 
nomial Ilpep °f 1- We might hope that in practice, the resolution of S 1 is 
relatively fast. 

3.4. A toric factorization algorithm. Proof of Theorem 1. By combining 
all previous results, we deduce the following 

Toric Factorization Algorithm (TFA). 

Input: f £ K[ti,ta] satisfying hypothesis {Hi) and (i^)- 
Output: The irreducible factorization / = q\ ■ ■ ■ q s of / over K. 

Step 0. Compute a regular fan E which refines S / . 

Step 1. Compute the set V of the irreducible rational factors of the exterior facet 
polynomials of /. 

Step 2. For i = 1, . . . , r and P £ Pi, compute the series <fip £ K[[x,]] with precision 
xf di using modular Newton iteration. 

Step 3. Compute the suitable powers of the 0p's in order to build the matrix A of 
Proposition 1. 

Step 4- Compute the reduced echelon basis of V(D). 

Step 5. Compute the rational factors of / by using Propositions 2 and 3. 

Theorem 1 follows immediately from the following 
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Proposition 4. The algorithm TFA is correct. It requires to factorize the exterior 
facet polynomials and to perform at most 0(Vo\(Nf) w ) arithmetic operations in K. 

Proof. The correctness of the algorithm is a consequence of Theorem 2 and Proposi- 
tions 1, 2, 3. The desingularization of the fan £ can be obtained by computing some 
Hirzebruch continued fractions (see [S]) and has a negligeable cost. We consider 
rational univariate factorization as a black-box of our algorithm. See for instance 
[15j . [21] and [4] for recent advances in that direction. Finally, the cost of steps 
2, 3, 4, 5 follows from Corollaries 1 and 2. □ 

4. Comparison with related results. Improvements 

In Subsection 4.1, we compare the algorithm TFA with the most related dense 
algorithms. In Subsection 4.2, we discuss the relation with the toric algorithm 
developed in [27) by the author. In particular, we obtain a sufficient criterion for 
using a smaller lifting precision. 

4.1. Comparison with dense algorithms. We compare our method with the 
lifting and recombination scheme proposed by Lecerf [TB] , [H] and Cheze-Lecerf [5] 
for dense polynomials and we discuss some possible improvements for each step of 
the algorithm TFA. 

About Step 0. Since C Pi Di — for all rays p, £ £ \ £/, we need not to compute 
all the fan £. Namely, we check easily that it's enough to compute the succesive 
rays pi+i £ £ of the rays pi £ £/. 

About Step 1. In most cases, the cost of the univariate factorization step dominates 
the complexity of the algorithm TFA. It might be interesting to avoid some of 
the facet factorization by chosing D with support strictly contained in \dX\. For 
instance, if / has bidegree (di,^), then X = P 1 x P 1 and we might hope to 
recover the decomposition of C from its restriction to D = {d\ + 1)P X since the 
corresponding restriction H°(X, Ox{C)) — > H°(X, Od{C)) is injective. This turns 
out to be the case : in [Tgj, G. Lecerf factorizes bidegree polynomials by using 
only one facet factorization with a sharp precision. In general, D has to obey to 
the vanishing cohomological properties used in the proof of Theorem 2, which are 
closely connected with the geometry of Nt. 

About Step 2. In the dense absolute case treated in [5], G. Lecerf and G. Cheze 
compute the analoguous series by introducing the Paterson-Stockmeyer evaluation 
scheme in the Newton iteration process. We can adapt such a method to our 
situation by replacing the input polynomial of Algorithm 1, [H] by a polynomial 
with degree bounded by Vol(iVA In such a way, the complexity 0(Vol(Nf) 2 ) of 
step 2 decreases to d(Vol(N f ) < - UJ + 1 ^ 2 ). 

About Step 3. In [TJ], G. Lecerf builds an analoguous linear system by using fast 
modular euclidean divison rather than by computing the 0p's powers. Both ap- 
proaches give rise to equivalent linear systems (see [5], Section 2.3), but the divison 
method permits to build the underlying matrix faster. We might hope that in the 
toric case, it is possible too to introduce an equivalent matrix that can be built 
using modular division. 
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About Step 4. In [18] and [19], the linear system resolution has complexity 0(d UJ+1 ) 
with d the total degree of /. It's easy to check that the sum of the lattice lengths of 
the exterior facets of Nf is bounded by d, with equality if and only if S/ is regular. 
It follows that Card('P) < d. Since C(Vol(iV/)) C 0{d 2 ), the reduced echelon basis 
computation is faster using the toric approach (much faster in the most case). 

About Step 5. Since we recover the factors of / by solving affine systems, step 5 has 
a relatively high cost in the algorithm TFA. If / is a dense polynomial, the task 
is much simpler and we can recover fastly the global factors of / from the local 
ones by using modular multiplications of (see |18j). or by using a partial fraction 
decomposition method (see [8], [1]). This permits a softly d 3 complexity for the 
factors computation, in general much faster than our approach. We might hope to 
adapt these methods to the toric case. 

4.2. About the lifting precision. In [19 , G. Lecerf gives an example in the 
dense case that shows that the precision D := 2div oc (/) is sharp in Theorem 2. 
On an other hand, the author obtain in [27j a toric factorization algorithm running 
with precision E := diVoo(/) +dX, but with exponential complexity in most cases. 
We explain here the relation with our algorithm and we give an explicit sufficient 
criterion for using the precision E < D. 

Let 7 = X)pe-p ^pIp £ V- We define the rational numbers 

hh) ■= Mpdeg(P), i=l,...,r 
Pev 

and we introduce the following convex subset of V 

r r 

A c = { 7 ey, ^efc.Tfc^ <£(e*,»*)*i(7) < 0, fc = l,2}, 

i=l z=l 

where (ei,e2) is the canonical basis. We have the following 

Proposition 5. The finite set V(E) n {0, 1} V n A c is a system of generators of 
the vector subspace V{D) C V(E). 



Proof. Since V(D) admits a basis with coordinates in {0,1}, it's enough to show 
that V{D) n {0, 1} V = V(E) n {0, l}l r l n A c . Let 7 e V(D) n {0, 1} V . So 7 
is restriction to D of a rational component C" £ Div(X) of C (Theorem 2). In 
particular, both divisors C and C—C are numerically effective, which is equivalent 
to that 

(15) < deg(C" • A) < deg(C • A), i = 0,...,r+l. 

For i = 1, . . . , r, we have equalities deg(C" • A) = ^(7) and deg(C • A) = h- O n 
an other hand, for any m £ Z 2 , we have 



r+1 r+1 

Y,(m,Vi)deg(C ■ A) =^(m,7 ?i )deg(C- A) = 0, 

i=0 t=0 
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since the divisor Yll=o ( m , Vi)Di — div(t m ) is principal. Letting m — e\ and using 
that (r) r +i,r)o) = (ei,e2), we deduce that 

r r 

deg(C'-D ) = - J>i,»fc)li(7), deg(C • D ) = - ]T(ei, m )l u 
i=l i=l 
and the same reasonning with m — e^ gives 

r r 

deg(C" • £> r +i) =-X)(ea»^(7). deg(C- - -£>2,^- 

i=l i=l 

Combined with (15), we deduce that 7 G Ac, giving an inclusion V^-D) fl{0, I} 73 C 

y(B)n{o,i} p nA c . 

Let us show the opposite inclusion. If 7 G V{E) (~l {0, 1} V n Ac, it lifts to 
some divisor C" 6 Div(Jf). By hypothesis, we have < 7 < 7c, giving obvious 
inequalities 

< deg(C • A) = h{l) <k = dcg(C • Di) 
for alH = 1, . . . , r. Since 7 G Ac, we deduce from the previous discussion that (15) 
holds. Since being nef is equivalent to being globally generated on a toric varietyit 
follows that both O x {C) and O x (C - C) are globally generated. By [23, proof 
of Theorem 2, this gives rise to supplementary vanishing cohomology properties 
which ensure that C can be chosen to be an absolute component of C. Since 7 is 
defined over K, such a component is rational by Lemma 8. Thus 7 is a {0, l}-linear 
combination of the 7/s, that is 7 G V(D) n {0, 1} V . □ 

Proposition 5 admits the following useful corollary, which is the toric version of 
[5], Proposition 4. 

Corollary 3. If all the exterior facets of Nf have an inward primitive normal 
vector with negative coordinates, then V(E) = V(D) if and only if each vector of 
the reduced echelon basis ofV(E) lies in {0,1} V . 

Proof. If the reduced echelon basis of V(E) lies in {0, we have ^(7) < h for 
all 7 in that basis. By hypothesis, we have (ek,i]i) < for alH = 1, . . . , r, k = 1, 2 
and it follows that 7 G V(E) n {0, \} v n A c . The equality V{E) = V(D) follows 
from Proposition 5. The other implication is trivial. □ 

Proposition 5 gives an efficient way to compute the reduced echelon basis of V(D) 
from the finite set V(E) n {0, 1} V . Roughly speaking, the underlying algorithm is 
that developed in [57] • It has the advantage to use a smaller lifting precision, but 
in return, it looks for "good" partitions of {(LI} 73 and can have an exponential 
complexity. We don't know what is the probability for that equality V(E) — V(D) 
holds. 

5. Conclusion 

We propose a new lifting and recombination algorithm for rational bivariate 
factorization that takes advantage of the geometry of the Newton polytope. For 
polynomials that are sparse enough, our complexity is competitive with that of 
the actual fastest algorithms developed for dense polynomials. We might hope 
to improve the complexity with a carefull application of the standard modular 
algorithms in the toric setting. 
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